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In typical plasma physics textbooks, Debye shielding, coUisional transport. Landau damping of 
Langmuir waves, and spontaneous emission of these waves are introduced in different chapters. This 
paper provides a compact unified introduction to these phenomena without appealing to fluid or 
kinetic models, but by using Newton's second law for a system of N electrons in a periodic box 
with a neutralizing ionic background. A rigorous equation is derived for the electrostatic potential. 
Its linearization and a first smoothing reveal this potential to be the sum of the shielded Coulomb 
potentials of the individual particles. Smoothing this sum yields the classical Vlasovian expression 
including initial conditions in Landau contour calculations of Langmuir wave growth or damping. 
The theory is extended to accommodate a correct description of trapping or chaos due to Langmuir 
waves. In the linear regime, the amplitude of such a wave is found to be ruled by Landau growth 
or damping and by spontaneous emission. Using the shielded potential, the coUisional diffusion 
coefficient is computed for the first time by a convergent expression including the correct calculation 
of deflections for all impact parameters. Shielding and coUisional transport are found to be two 
related aspects of the repulsive deflections of electrons. 

I. INTRODUCTION 

From the outset, inspired by gas physics, plasma physicists derived kinetic equations to describe microscopic aspects 
of their physics, in particular the Vlasov equation. This trend has been the dominant one till nowadays. However, 
for plasmas where transport due to short range interactions ( "collisions" ) is weak, it is possible to work directly with 
TV-body classical mechanics. As will be recalled in section |Vl this approach led to a new description of wave-particle 
interaction making it more intuitive and unifying particle and wave evolutions, as well as collective and finite-iV 
physics. 

This paper develops this approach further by deriving a rigorous fundamental equation for the electrostatic potential 
(section HIT)) , which is extended to enable a correct description of trapping or chaos due to Langmuir waves fsectionlV]). 
It brings further unifications: Debye shielding (or screening) and Landau theory of linear waves (section |IV|) . Landau 
growth or damping and spontaneous emission (section |V|, Debye shielding and coUisional transport (section IVip . 
These results come with a considerable simplification of the mathematical framework with respect to textbooks (see 
for instance Refs. [1H4|), and with new insights into microscopic plasma physics. In particular it solves the following 
issue: How can each particle be shielded by all other ones, while all the plasma particles are in uninterrupted motion? 
This turns out to be a mere consequence of the almost independent deflections of particles due to the Coulomb force. 
Furthermore, the new path to Landau damping goes first through Debye shielding, a totally unexpected fact, since 
classical textbooks present these concepts in different and unrelated chapters. 

As to coUisional transport, this paper does not simply recover the classical result presented in textbooks, but it 
provides the first calculation of this transport without any ad hoc cut-off. This is made possible by a mechanical 
calculation of the deflection of a particle which covers the scales corresponding to the Debye length, to the inter- 
particle distance, and to the distance of minimum approach of two electrons in a Rutherford collision. Till now none 
of the papers about coUisional transport, and no textbook, was able to correctly describe the scale of the inter-particle 
distance. 

The derivations are elementary (we refrain from invoking the Liouville equation or the BBGKY hierarchy). After 
section IIIIl the successive approximations of the fundamental equation (jlSp enable to get rid of the heavy notations 
related to its generality. These approximations involve linearization in sections IIIIl IIV Al and IIV CI but section lIVBi 
the beginning of section |Vl and section |VI] do not appeal to linearization. For the sake of brevity, this paper is more 
oriented toward concepts and intuitive physics, than toward ready-to-use formulas to be presented elsewhere. 

II. CLAIMS 

Here are the main results of this paper: 

1. By using the Fourier and Laplace transforms, a rigorous equation (Eq. ([T5|) ) is derived for the electrostatic 
potential of an infinite plasma made up of the periodic replication of N electrons coupled by Coulomb forces in a 
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volume with a neutralizing ionic background (One Component Plasma (OCP) model 0-01 )• This equation is of 
the type Eip = source term, where f is a linear operator, acting on the infinite dimensional array Lp whose components 
are all the Fourier-Laplace components of the potential. 

2. Considering that particle trajectories stay close to ballistic, and replacing inside S the discrete sums over 
particles by integrals over a smooth distribution function, the potential turns out to be the sum of the shielded 
Coulomb potentials of the individual particles (Eq. ()24p ). Such potentials were first computed by a kinetic approach 
in section II. A of Ref. fi| and later on in [1, 

3. Replacing the discrete sums over particles of their shielded potentials by integrals over a smooth distribution 
function, yields the classical Vlasovian expression (Eqs (|20t and (|33|) ) including initial conditions in Landau contour 
calculations of Langmuir wave growth or damping, obtained by linearizing Vlasov equation and using Fourier-Laplace 
transform, as described in many textbooks (see for instance Refs. [J, |3, Q)- Therefore in these calculations, the 
electrostatic potential turns out to be the smoothed version of the actual shielded potential in the plasma. 

4. In the spirit of Refs. [TTI - [T3 |. in order to accommodate a correct description of trapping or chaos due to Langmuir 
waves, the set of particles is split into bulk and tail, where the bulk is the set of particles which cannot resonate with 
Langmuir waves. Repeating the analysis leading to Eq. ()15|) for the bulk particles, the same equation is recovered 
with an additional source term due to the tail particles (Eq. p6p ). 

5. Using the fact that the number of tail particles is small with respect to the bulk one, and a technique introduced 
in Refs. [isl. [l^. an amplitude equation is derived for any Fourier component of the potential where tail particles 
provide a source term (Eq. ([iP]) ). 

6. This equation, together with the equation of motion of the tail particles, enables to show that, in the linear 
regime, the amplitude of a Langmuir wave is ruled by Landau growth or dam ping, a nd by spontaneous emission 
(Eq. (I42p ). a generalization to 3 dimensions of a one-dimensional result of Refs. |12l. [l5|. 

7. Using the shielded potential, the collisional diffusion coefficient 13 of a given particle is computed by a convergent 
expression including the particle deflections for all impact parameters. These deflections are computed by flrst order 
perturbation theory in the total electric field, except for those due to close encounters. The contribution to D of the 



former ones is matched with that of the latter ones provided by Ref. 16 | . The detailed matching procedure includes 
the scale of the inter-particle distance, and is reminiscent of that in 17|, without invoking the cancellation of three 
infinite integrals. D has the same expression as that in Ref. (l6j . except for the Coulomb logarithm which is modified 
by a velocity dependent quantity of order 1 . 

8. Picard iteration technique is applied to the equation of motion of particle P due to the Coulomb forces of all 
other ones. It reveals that a part of the effect on particle P of another particle P' is mediated by all other particles 
(Eq. ([31]) ). Indeed particle P' modifies the position of all other particles, implying that the action of the latter ones 
on particle P is modified by particle P' . 

9. This calculation yields the following interpretation of shielding. At t = consider a set of randomly distributed 
particles, and especially particle P. At a later time i, the latter has deflected all particles which made a closest 
approach to it with an impact parameter b ^ Vt^t where wth is the thermal velocity. This part of their global 
deflection due to particle P reduces the number of particles inside the sphere S(t) of radius Vtht about it. Therefore 
the effective charge of particle P as seen out of S{t) is reduced: the charge of particle P is shielded due to these 
deflections. This shielding effect increases with t, and thus with the distance to particle P. As a result, the typical 
time-scale for shielding to set in when starting from random particle positions is the time for a thermal particle to 
cross a Debye sphere, i.e. LiJ~^^ where cjp is the plasma frequency. Furthermore, shielding is a cooperative dynamical 
process: it results from the accumulation of almost independent repulsive deflections with the same qualitative impact 
on the effective electric field of particle P (if point-like ions were present, the attractive deflection of charges with 
opposite signs would have the same effect). It is a cooperative effect, but not a collective one. So, shielding and 
collisional transport are two aspects of the same two-body repulsive process. 



III. FUNDAMENTAL EQUATION FOR THE POTENTIAL 

This paper deals with the One Component Plasma model 0-01 ) which considers the plasma as infinite with spatial 
periodicity L in three orthogonal directions with coordinates {x, y, z), and as made up of N electrons in each elementary 
cube with volume L^. Ions are present only as a uniform neutralizing background, which enables periodic boundary 
conditions. This choice is made to simplify the analysis which focuses on (p{r), the potential created by the N particles 
at any point where there is no particle. The discrete Fourier transform of ip(r), readily obtained from the Poisson 
equation, is given by ifiO) = and, for m 7^ 0, by 
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where — e is the electron charge, Eq is the vacuum permittivity, is the position of particle j, S is the set of integers 
going from 1 to N, (p{m) — J (^(r) exp(— ikf„ • r)d'^r, m ~ {mj-^my^niz) is a vector with three integer components, 



km = x"™' ^"^^ ^™ ~ l|l^m||- Reciprocally, 



m 

where X^m means the sum over all components of m running from — oo to +cx). 
The dynamics of particle I is determined by Newton's equation 

n = — V^,(n), (3) 

where mg is the electron mass, ipi is the electrostatic potential acting on particle i.e. the one created by all other 
particles and by the background charge. Its Fourier transform is given by Eq. ([T]) with the supplementary condition 
j ^ I. Let 

r["^ = r,o + ^it (4) 

be a ballistic approximation of the motion of particle /. In the following, we consider two instances of this approx- 
imation: the one where rjo and v; are respectively the initial position and velocity of particle Z, and the one where 
they are slightly shifted from these values by low amplitude Langmuir waves. Let 5vi = r; — r[°\ so that Eq. ^ is 
equivalent to 

5vi = V' kn ^/(n)exp[ikn • (r|°^ + Sri)]. (5) 



and 



Sections IIVI and fVl consider cases where the Jr^'s are small in some sense. To this end, we split (fiim) as 

<^;(m) = (f>i{m) + Aif}i{m.) (6) 

where 
with 

S4)j{m.) = ^exp(-ikm • ^ ikm • Svj), (8) 

A<^,(m) = ^ V exp(-ikm.rf )i?,(m), (9) 

with 

i?,j(m) = exp(-ikm • ^r^) - 1 + ik^ • Srj, (10) 

which is of order two in Srj . 

We now introduce the time Laplace-transform which transforms a function f{t) into /(w) = f(t) exp(icjt)dt 
(with Lu complex). The Laplace transform of Eq. ([5]) is 

w^5f;(cj) = - ^3^ ^knexp(ikn • r;o) "ifiiipi -,11,10 + ujnA) + iujSri{0) - Sri{0). (11) 

n 

where carets indicate the Laplace transformed versions of the quantities in Eq. ([5]), uj^.i = k„ • v; comes from the time 
dependence of r[°^ in the exponent of Eq. ([5]), and the operator acting on a function g{in,uj) is defined by 

*,(5;n,.)=ff(n,-)*r,(n,-), (12) 
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where • stands for the frequencies, * is the convolution product, and T;(n, w) is the Laplace transform of exp(ik„ • 5ri). 
The Laplace transform of Eqs (p|)- pll)) yields 

fc^(p;(m, = fc^0|°°^(m, + — exp(-ikm ■ Tjo) [km ■ (5fi(t^ - Wmj) + i^j(m,a; - Wm.j)], (13) 

where Rj{m,uj) is the Laplace transform of Rj, and (f)[^^\m,uj) is the Laplace transform of (^((m) computed from 
Eqs (O and (|S]) by setting Srj — for all j's in the latter. Substituting the Srj^s with their expression Eq. pT|) yields 

fcL^Km,c.)--3^^km.k„ 5: ' " + ""'^r ""'■^•^ exp[i(k„ - km) . r,o] 

= k'i,(j)f\ra,uj) - — V exp(-iki„ • rjo)-Rj(a; - Wmj), (14) 

where (m, w) is the Laplace transform of 4>i (m) computed from Eqs (O and (H)) by setting now Svj ~ Srj {0)+drj {0)t 
for all j's in the latter. 

Summing Eq. (jl4p over I — 1, ...N and dividing by — 1, yields 

fcmV(m, uj) - -3 2^ km • k„ 2^ -2 exp[i(k„ - k^) • r^o] 

= kl^^^"\m,uj) - — V'exp(-ikm • rjo)^j(w - Wmj), (15) 

where 4>^'^\m,uj) is (^["•^ (m, w) complemented by the missing l-th term. Equation ((T5|) is the fundamental equation 
of this paper, and is a rigorous consequence of Eqs ([T]) and ([3]): no approximation was made. This fundamental 
equation is of the type Sip = source term, where f is a linear operator, acting on the infinite dimensional array whose 
components are all the (^(m, aj)'s. 



IV. DEBYE SHIELDING 
A. Shielded Coulomb potential 

We now specialize Eq. (|15l) by considering the lowest order contribution of the i^r^'s, which makes the Rj's vanish 
(Approximation 1). We further consider (p to be small, and the (5r;'s to be of the order of (p (Approximation 2). These 
two approximations reduce ^'j ((^; n, w) to ip{n,uj). This corresponds to substituting the true dynamics in Eq. ([3]) with 
an approximate one ruled by 

Sri^—\/Mr\°^+Sri), (16) 

771 e 

where 0;(r) — J2j£S-j^i ^^'ji^) inverse Fourier transform of Eq. so that 

e eSrj ■ (r — rf^) 

hm S^,{r) = '-^ -f^, (17) 

47reo||r-rf II 4^eo||r - f 

The jth component of the approximate electric field acting over particle I turns out to be due to a particle located 
at rj°' instead of r^, and is made up of a Coulombian part and of a dipolar part with a dipole moment —eSrj. The 

cross-over between the two contributions occurs for ||ri — r^'^'^ll on the order of ||(5rj||, i.e. when the distance between 
particle / and the ballistic particle j is equal to the distance between the latter and the true particle j. For larger 
values of ||ri — rj^^H, the dipolar component is subdominant. For smaller ones, it is dominant, but with a direction 

which is a priori random with respect to the Coulombian one ((r; — rj'^'') is almost independent from Svj). Since the 
\\Srj II 's are assumed small, the latter case should be rare since it corresponds to a very close encounter between particle 
I and the ballistic particle j. As a result the approximate electric field stays dominantly of Coulombian nature, but 
with a small mismatch of the charge positions with respect to the actual ones. 
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We introduce a smooth function /(r, v), the smoothed velocity distribution function at < = such that the distri- 
bution 



les 



»f{r,v)d^rd^v + W{»), (18) 



where the distribution W yields a negligible contribution when applied to space dependent function which evolve 
slowly on the scale of the inter-particle distance; there the spatial integration is performed over the elementary cube 
with volume L'^, and the velocity integration is over all velocities. 

Replacing the discrete sums over particles by integrals over the smooth distribution function /(r, v) (Approximation 
3), Eq. ([T5|) becomes 

klHm, 0.) = fc^0(o) (m, ^) + ^ k„ ■ k„ / ^K^ + (kn-k )-v) ^-^^ _ ^ ^3^^ ^^^^ 

Li m^eo ^ J [Lu Km • vj 

where $ is the smoothed version of (p resulting from Approximations 1 to 3, and / is the spatial Fourier transform 
of /. We further assume the initial distribution / to be a spatially uniform distribution function /o(v) plus a small 
perturbation of the order of $ (in agreement with Approximation 2). Then operator £ becomes diagonal with respect 
to both m and ut (a complex quantity), and linearizing Eq. for <t yields 

e(m,tj)<l(m,w) = <^(")(m,a;), (20) 

where 

e(m, ..) = 1 - ( y d^v. (21) 

L'^meeo J (w - km • v)^ 

This shows that the smoothed self-consistent potential $ is determined by the response function e{m,uj). It is the 
classical plasma dielectric function. A first check of this can be obtained for a cold plasma: then e(m, w) = 1 — Wp/w^, 
where Wp = [(e^n)/(meeo)]^/^ is the plasma frequency {n = N/L^ is the plasma density). The classical expression 
involving the gradient of /o in v is obtained by a mere integration by parts. 

The contribution of particle j to ^("-'(m) is (5^j''''(m) = — ^J^^ exp[— ikm • (fjo + Vji)]. Its Laplace transform is 

.7(0), ^ ie exp(-ikm • r^o) 

6(^)'{m,uj) = 22) 

The corresponding part of $(m,w) is 5$j(m,w) = (50^°''(m,w)/e(m, cj). This turns out to be the shielded potential 



of particle j [SHlOi . By inverse Fourier-Laplace transform, after some transient whose duration is estimated at the 
end of section IIVBI the potential due to particle j becomes the shielded Coulomb potential 

<5$,(r)=<5<f(r-r,o-v,t,v,), (23) 

where 

i^eo ^fc^e(m,km-v) 

Therefore, after this transient, the dominant contribution to the full potential in the plasma turns out to be the 
sum of the shielded Coulomb potentials of individual particles located at their ballistic positions. Let Ad = 
[{eoksT) / (ne^)]-^^^ ~ [{k-BT) /rriel^^^ujp^ be the Debye length, where fce is the Boltzmann constant, T is the tem- 
perature, and n is the density. The wavenumbers resolving scale ||r|| are such that fcm||i"|| ^ 1. If ||i"|| <C Ad, 
the corresponding wavenumbers are such that /cmAD 3> 1. Therefore there is no shielding for ||r|| <C Ad, since 
e(m, km • v) - 1 ~ (fcmAD)"^. 

In the following, Eq. (^5)) is used by substituting (5$(r — r^o — Vjt, Vj) with 5$(r — r^, Vj): the shielded potential 
of particle j is computed by taking into account its actual position, since it is the original Coulomb one close to rj. 
The error made for r — Vj of the order of Ad is small as long as the mismatch of from the ballistic orbit is much 
smaller than Ad • As was done for the bare potential of Eq. ([Ij , the field acting on a given particle I is obtained by 
removing its own divergent contribution 6^i from $. 
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B. Mediated interactions imply Debye shielding 



In the above introduction of Debye shielding, using the Laplace transform of the particle positions does not provide 
an intuitive picture of this effect. We now show that such a picture can be obtained by turning back to the mechanical 
description of microscopic dynamics with the full Coulomb potential of Eq. ([T]). In order to compute the dynamics, 
we use Picard iteration technique. From Eq. rl"-*, the nth iterate of r^, is computed by 

r^r^ = ^Wr'\r\-~'^ (25) 

where (p[" is computed by the inverse Fourier transform of Eq. (IlJ with the Vj's substituted with the r^" ^■''s. The 
iteration starts with the ballistic approximation of the dynamics defined by Eq. Q , and the actual orbit of Eq. ([3]) 
corresponds to n oo. Let (5r|"'' = r|"'' — rf^^ be the mismatch of the position of particle I with respect to the 
ballistic one at the nth iterate. It is convenient to write Eq. (|25|) as 

S^\-^^ Y: 5^. (26) 



with 



and 



<5t[;'=.ac(r|"-^'-r("-^)), (27) 



ac(r) = V^(pc(r), (28) 



with ^i^C; the L-periodized, OCP bare Coulomb potential, given by 

5^c{v) - exp(ik™ • r), (29) 

such that (p(r) = Yj^es ^Vc{^ - ^j)- Then 6v^j^^ = J^jes-.^^i ^c(r|"'' - rf^) and for n > 2 

Sr\"^ - Sr['^ + Y Vac(rJ"' - rf ) • {6r\"-^^ - Jrf-^^) + 0{a^), (30) 

where a is the order of magnitude of the total Coulombian acceleration. Equation pop may be written 

f|"^= E P?g^+<-'V2Vac(r|°)-rf).<-^V0(a^), (31) 

where Mj^ = Vac(r['^'' — v^f^) ■ J2i£S-iiii j^^^ii^ ~ ^^'ji ^'') modification of the bare Coulomb acceleration 

of particle j on particle / due to the following phenomenon: particle j modifies the position of all other particles, 
which implies the action of the latter ones on particle I is modified by particle j. Therefore m/J' is the acceleration 
of particle I due to particle j mediated by all other particles. The last term in the bracket in Eq. (|3T|) . where 
(5r|"~"'^-' = /o /o Sr\j~^\t")dt"dt' , accounts for the fact that both particles j and I are shifted with respect to their 
ballistic positions. 

Since the shielded potential of section HV Al was found by first order perturbation theory, it is felt in the acceleration 
of particles computed to second order. This acceleration is provided by Eq. pip for n = 2. Therefore its term in 
brackets is the shielded acceleration of particle I due to particle j. As a result, though the summation runs over all 
particles, its effective part is only due to particles j typically inside the Debye sphere about particle I. Starting from 
the third iterate of Picard scheme, the effective part of the summation in Eq. (I5T|) ranges inside this Debye sphere, 

since the 5r["'~^'''s are then computed with a shielded acceleration. This justifies the use of the shielded potential to 
compute collisional transport in section IVTl and is the basis of the intuitive interpretation of Debye shielding in claim 
9. This interpretation shows that, when starting form random particle positions, the typical time-scale for shielding 
to set in is the time for a thermal particle to cross a Debye sphere, i.e. uj~^, which sets the duration of the transients 
occurring in the inverse Laplace transform leading to Eq. (|23p . This order of magnitude is correct for a plasma close 
to equilibrium. 
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C. Debye shielding and Landau damping 

We now apply the smoothing using distribution function / to 0^"^(m, w) too in Eq. ((20|) . As a result of Eqs ([TllSl), 
this yields 

e 



l>(")(m) = -—^ 1 1 exp[-ik„ • (r + vi)]/(r,v)d3rd'V (32) 



whose Laplace transform is 



Eofcm J - km ■ V 



(33) 



which shows this second smoothing makes Eq. to become the expression allowing Landau contour calculations, 
as stated in Claim 3. 

We point out that, in this paper, the smoothed velocity distribution is introduced after particle dynamics has been 
taken into account, and not before, as occurs when kinetic equations are used. This avoids addressing the issues of 
the exact definition of the smoothed distribution for a given realization of the plasma, and of the uncertainty as to 
the way the smoothed dynamics departs from the actual iV-body one. 

V. WAVE-PARTICLE DYNAMICS 

In section llVCi the existence of Langmuir waves is asserted by connection with Landau theory. To describe 
Langmuir waves with discrete particles, we consider that the rjo's are random, and we allow for non zero 5rj(0)'s and 
(5rj(0)'s for the ^r^-'s in Eq. (|8]). Therefore, in the formulas of section HVl the r^o and Vj's are slightly shifted with 
respect to the initial rj(0)'s and rj(0)'s due to Langmuir waves. 

Up to this point we described Langmuir waves by a linear theory. We now generalize the analysis of section IIIII 
to afford the description of nonlinear effects in wave-particle dynamics. Indeed, resonant particles may experience 
trapping or chaotic dynamics, which imply km • Jrj's of the order of 27r or larger for wave km's. To describe such a 
dynamics, it is not appropriate to expand (j) as was done in Eqs ([TllS]) for such particles. However, this expansion may 
still be justified for non resonant particles over times where trapping and chaos show up for resonant ones. In order 
to keep the capability to describe the latter effects, we now split the set of N particles into bulk and tail, in the spirit 
of Refs. |ll|-[l5|. The bulk is defined as the set of particles which are not resonant with Langmuir waves. We then 
perform the analysis of section Hill for the iVbuik particles, while keeping the exact contribution of the iVtaii particles 
to the electrostatic potential. To this end, we number the tail particles from 1 to iVtaii, the bulk ones from iVtaii + 1 
to = A^buik + A'taii, and we call these respective sets of integer S'taii and 5'buik- For I G ^buik, we now substitute 
Eq. dH) with 

^Km) = ^';""^~V (m)+ J2 ^'^^■("^)' (34) 

where 



jeSbuiky'^i 



t/(m) = ^t^r^^ — TT exp(-ikm • r^). (35) 

£0fcm(^bulk - 1) ^^^^^ 



In the r.h.s. of Eq. ((M)) . the first term vanishes if A^taii = 0. We now perform the calculation of section IIIII on 
substituting the previous summations with index running from 1 to A'^ by ones where the index runs over S'buik, while 
keeping the exclusion oi j — I where indicated. The previous division by A^ — 1 preceding Eq. (jlSp is now a division 
by A^buik - 1- This yields 

kL^M E^^-^r^Y: ' " + ""^'^"^ ^^P[i(kn - km) ■ r,o] 

n j£S ^ '■' ' 

= /c^0("'(m, w) - — V'exp(-ikm • rjo)i?j(tj - Wmj) + A:m?7(m,a;), (36) 

where L/(m,cj) is the Laplace transform of ?7(m). Then Eq. ((20)) becomes 

= 0(°)(m,a;) + (37) 
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Let $(m, be the inverse Laplace transform of $(m,a;), ^buik(m, w) be the solution of Eq. ([20)) computed for the 
bulk particles, and j>hnik(m, t) be its inverse Laplace transform. We now derive an amplitude equation for $(m,i) 
in a way similar to Refs. [3, Let Wm be such that e(m, Wm) = 0; because of the definition of the bulk, this 
frequency is real. Then $buik(m,0 — ^exp(— iwmOi where A is a constant, and 

^(o)(m,^) = ^^, (38) 

according to Eq. (1^ . 

Let g{m,t) = $(m, i)/$buik(m, i). Therefore <l>(m,a;) = Ag{uj-uJm), which together with Eq. ([37]) and dSS]) yields 



Ae{m,uJm+uj')[g{m,Lu') - —] = U{m,uj^ + uj'), (39) 



where uj' = lo — LOm- If -^taii -^buik, 5(m, t) is a slowly evolving amplitude, and the support of g(m, w) is narrow 
about ze 
order. S 
l>(m, t) 



about zero. This justifies Taylor-expanding e(m, Wm+i^') about = in Eq. ([M)) . which yields Ms^fuii^;' to lowest 
order. Setting this in Eq. p9p and performing the inverse Laplace transform finally yields an amplitude equation for 



d'^{m,t) ~ ieiVbuik / -i ^ fAn\ 

gofc^(A^buik - 1) ^ a: j-gg^^„ 

The self-consistent dynamics of the potential and of the tail particles is ruled by this equation and by the equation of 
motion of these particles 

^k„ l-,(n)exp(ik„T,). (41) 



' Lh 

n 

These two sets of equations generalize to 3 dimensions the self-consistent dynamics defined in Refs. [ill [13 • ^'^^ 
the sake of brevity, we do not develop here the full generalization of the analysis in these papers; it is lengthy, but 
straightforward. However, since this analysis unifies spontaneous emission and Landau growth and damping, we give 
the result ruling the evolution of the amplitude of a Langmuir wave provided by perturbation calculation where the 
right hand sides of Eqs (|40ll4ip are considered as small of order one. This is natural for Eq. (|40p since iVtaii ^ ^buik, 
and for Eq. (pij) if the Langmuir waves have a low amplitude. Let J(m,i) — ($(m, t)$(— m, t)), where the average 
is over the random initial positions of the tail particles (thus their distribution is spatially uniform). Then a second 
order calculation in $ yields 

^^^^^^ = 27™l J(m, t) + 5„,,po„t, (42) 
where 7mL is the Landau growth or damping rate given by 

7mL=am— 1 — [-, — ;m) (43) 
av km 



with 



2 

7re 



m, 



and /rcd('y;iii) is the reduced smoothed distribution function /i.od(w;ni) — ///(wkm + vj^) d'^v± where km is the 
unit vector along km and vx is the component of the velocity perpendicular to k^; 5'mspont is given by 

spoilt — 2/ Jrcdy j 7' V '-^/ 



where n = N jl? is the plasma density. S'm spont corresponds to the spontaneous emission of waves by particles and 
induces an exponential relaxation of the waves to the thermal level in the case of Landau damping. The second order 
calculation for the particles yields the diffusion and friction coefficients of the Fokker-Planck equation ruling the tail 
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dynamics. This equation corresponds to the classical quasilinear result, plus a dynamical friction term mirroring the 
spontaneous emission of waves by particles, as found in the one-dimensional case in Refs. [l^.[l5|. 

Finally, we point out that Eqs (|40II41|) enable to tackle the linear theory of Langmuir waves by considering a single 
mechanical system [l5j . This is naturally obtained by considering the unperturbed plasma as made up of a series of 
monokinetic parallel beams, and each beam as an array of equidistributcd particles. This approach finds eigenmodes 
corresponding to the Landau instability [3] , if the distribution function has a positive slope. If the slope is negative, 
it recovers the van Kampen modes 'l9l|, as a series of beam modes providing Landau damping as a result of the 
phase mixing of the latter. This phase mixing turns out to be also important to recover a correct growth of the wave 
amplitude in the unstable case (section 3.8.3 of reference We point out that in Refs. [H, [HI the equivalent of 

Eqs (|40l4ip was obtained without using any smoothing, but by a mechanical reduction of degrees of freedom, which 
made even stronger the idea of considering a single mechanical system. 



VI. DEBYE SHIELDING AND COLLISIONAL TRANSPORT 



We now focus on the case where the particles have random initial positions, i.e. where the plasma has a uniform 
density, and for simplicity we consider the plasma to be in thermal equilibrium. Then the dynamics of particles has no 
collective aspect, but is ruled by the cumulative effect of two-body deflections. More specifically, we choose random 
r/o's, and vanishing (5r/(0)'s and (5r/(0)'s, and we assume that at t = all particles are in the same cube with volume 

, and we consider the limit L/Ad — > oo. In contrast, each particle has a well defined velocity, in such a way that the 
overall initial smoothed velocity distribution is close to some given function. We focus on particle I which is assumed 
to be close to the center of the cube. In this section, we approximate the true dynamics by that due to the shielded 
Coulombian interactions, i.e. we write 

5yi^ H - ''i'Vj)' (46) 

with 

a(r,v) = — V(5$(r,v), (47) 

me 

where (5$(r,v) is given by Eq. ([M)) . 

We compute particle / deflection in a sequence of steps. First, we use flrst order perturbation theory in which 
shows the total deflection to be the sum of the individual deflections due to all other particles. For an impact 
parameter much smaller than Ad, the deflection due to a particle turns out to be the perturbative value of the 
Rutherford deflection due to this particle if it were alone. Second, for a close encounter with particle n, we show that 
the deflection of particle I is exactly the one it would undergo if the other N — 2 particles were absent. Third, the 
deflection for an impact parameter of order Ad is given by the Rutherford expression multiplied by some function of 
the impact parameter reflecting shielding. These three steps yield an analytical expression for deflection whatever the 
impact parameter. 

We first compute 6ri by flrst order perturbation theory in taking the ballistic motion defined by Eq. as 
zeroth order approximation. This yields 

6rn{t)^ J2 SrijiiO,t), (48) 

where 

Srijiih,t2) = /*' a[rj'')(i') - rf\t'),^,]dt' . (49) 

It is convenient to write 

r['\t') - rf \t') = b,, + Avy (i' - U,), (50) 

where Uj is the time of closest approach of the two ballistic orbits, and hij is the vector joining particle j to particle 
I at this time. Then bij — \\hij\\ is the impact parameter of these two orbits when singled out. The initial random 
positions of the particles translate into random values of hij and of tij. The typical duration of the deflection of 
particle I given by Eq. (j49l) is Atij = bij/Avij where Awjj = ||Av;j||, but a certain number, say a, of Atjj's are 
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necessary for the deflection to be mostly completed. For a given bij and for t ^ Atij , the deflection of particle I given 
by Eq. (|49p is maximum if is in the interval [aAtij,t — aAtij]. We notice that Atij is about the inverse of the 
plasma frequency for bij ^ Ad and Ai;;^ on the order of the thermal velocity. 

For the sake of brevity, we compute here just the trace of the diffusion tensor for the particle velocities. To this 
end, we perform an average over all the r/o's to get 

{Srlit))^ J2 KiW), (51) 

taking into account Eq. (|24)) . and the fact that the initial positions are independently random, as well as the — r^'s 
for i ^ j. Therefore, though being due to the simultaneous scattering of particle I with the many particles inside 
its Debye sphere, (Srfi{t)) turns out to be the sum of individual two-body deflections for 6;j's such that first order 
perturbation theory is correct. Hence the contribution to {Sirfi {t)} of particles with given bij and Avij can be computed 
as if it would result from successive two-body collisions, as was done in Ref. flG*] and in many textbooks. 

For an impact parameter much smaller than Ad, the main contribution of a[r|°'' (t') — r'"^^ (t'), Vj] to the deflection of 
particle I comes from times t' where ||r|"''(t') — rj"''(t')|| <^ Ad- Therefore a(r, v) takes on its bare Coulombian value, 
and ((5rfi(i)) is a first order approximation of the effect on particle I of a Rutherford collision with particle j. Comparing 

2 

this approximate value with the exact one shows the perturbative calculation to be correct for bij 3> A^a = sqAv'^- ' 
the distance of minimum approach of two electrons in a Rutherford collision, as given by energy conservation. 

Second, we consider the case of the close approach of particle n to particle I, i.e. bin ~ Ama- We write the acceleration 
of particle I as 

r;=a(r;-r„)+ ^ a(r/-r-,). (52) 

For particle n, we write the same equation by exchanging indices I and n. Since the two particles are at distances 
much smaller than the inter-particle distance d — n^^^^ = L/N^^^, the accelerations they get from all other particles 
are almost the same. Therefore, when subtracting the two rigorous equations of motion, the two summations over j 
almost cancel, yielding 

d2(r, - r„) 

— = 2a(r, - r„), (53) 

which is the equation describing the Rutherford collision of these two particles in their center of mass frame, in the 
absence of all other particles (at such distances the shielded potential is the bare Coulomb one). Since bin ^ d, At in 
is much smaller than the Atij's of the other particles. Therefore the latter produce a negligible deflection of the center 
of mass during the Rutherford two-body collision, and the deflection of particle I during this collision is exactly that 
of a Rutherford two-body collision. The contribution of such collisions to {Srf{t)} was calculated in Ref. [l6l |. 

Now, since the deflection of particle I due to particle j as computed by the above perturbation theory is an 
approximation of the Rutherford deflection for the same impact parameter, we may approximate the perturbative 
deflection one by the full Rutherford one, which provides an obvious matching of the theories for bij A^a and for 
Ad ^ bij ^ Ama: we may use the estimate of [iB] in the whole domain bij <C Ad- 

Third, we deal with impact parameters of the order of Ad. Then the deflection due to particle j must be computed 
with Eq. (j49p . For the sake of simplicity, we make the calculation for the case where Vj is small, which makes 
<5$(r,v) ~ 5$(r,0) which is the Yukawa potential 5$Y(r) = -4:^;^ exp(-M) (Eq. (18) of Ref. 0). The first 
order correction in km • Vj to this approximation is a dipolar potential with an electric dipole moment proportional 
to Vj . Since a Maxwellian distribution is symmetrical in v, these individual dipolar contributions cancel globally. As 
a result, the first relevant correction to the Yukawa potential is of second order in km ■ Vj. This should make the 
Yukawa approximation relevant for a large part of the bulk of the Maxwellian distribution. 

In the small deflection limit, a calculation using the fact that the force derives from a central potential shows the 
full deflection of particle I due to particle j is provided by 

e2 f+°° 1 1 r(t) 

5vijii-oo,+oG) =. b,, / [—- + ]exp[— ^]dt, (54) 

47rmeeo J^oo (*) AD?"^(i) Ad 

where r{t) — [bf- + Avf-t^Y^"^ and hij was defined with Eq. ([50]) . Defining 9 = aics\a[Avijt / r{t)], this equation 
becomes 



2e2 h{bi£) 



(5rzji(-oo, +00) = "2 — \iij, (55) 

47rmeeoA'yi,- bf, ' 
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where 

hib) = r^\cos{e) + ^]exp[- \ ] de<[l + ^]exp[-^]. (56) 
Jo -^D ADCOs(t^) 2Ad Ad 

During time t ^ ^tij, a volume 2T:AvijtbijSbij of particles with velocity Vj and impact parameters between bij and 
bij + dbij produce the deflection of particle I given by Eq. ([55]) . and a contribution scaling like —^p^Sbij to {Srfi(t)). 
Let &,„iii be such that Ad 3> 6min 3> Ama- The contribution of all impact parameters between 6min and some &max 
is thus scaling like the integral /jl'"'''' h?{b)/b db. Since for b small /i(0) ~ 1, if 6max ^ Ad, this is the non-shielded 
contribution of orbits relevant to the above perturbative calculation. Since, by approximating it by the Rutherford-like 
result of Ref. [l6| , this contribution matches that for impact parameters on the order of Ama, the contribution of all 

impact parameters between A,na and some 6max small with respect to Ad is thus scaling like the integral /^°""' \/b db 
as was computed in Ref. [l6l |. The matching of this result for 6 ^ Ad is simply accomplished by setting a factor 
K^il)) in the integrand which makes the integral converge for b — t- oo. Taking this limit, one finds that the Coulomb 
logarithm ln(AD/A„ia) of the second Eq. (14) of Ref. 116] becomes ln(AD/Ama) + C where C is of order unity. If the 
full dependence of the shielding on Vj were taken into account, the modification of the Coulomb logarithm would be 
velocity dependent. 

For the sake of brevity, we do not develop here the calculation of the dynamical friction which will be presented 
elsewhere. This calculation requires using second order perturbation theory, but follows the same lines as those for 
the diffusion coefficient. 

For an inhomogeneous plasma, the acceleration of particle I may be split into a homogeneous and a wave part. 
When using the linearized versions of Eq. (|15p and subsequent ones in section HVl we can split all (/)'s and $'s into 
a homogeneous part and an independent inhomogeneous one. Therefore the diffusion coefficient and the dynamical 
friction estimated by perturbative calculation of the dynamics up to second order are the sum of the collisional 
contribution and of a contribution due to waves, the latter as calculated for instance in Refs. [13, [3 ■ For the sake 
of brevity, we defer this point to later publication. 



VII. CONCLUSION 



In this paper, Debye shielding, collisional transport, Landau damping of Langmuir waves, and spontaneous emission 
of these waves were introduced by a mechanical approach. This approach brings unification and simplification in basic 
microscopic plasma physics, and may be useful for pedagogical purposes. One might think about trying to apply the 
above approach to plasmas with more species, or with a magnetic field, or where particles experience trapping and 
chaotic dynamics. The first generalization sounds rather trivial, and the third one is under way, at least in one 
dimension (see a pedestrian introduction in [20'| and more specific results in pll [2^). 

As in many textbooks, linearization was applied in sections IIIIl IIV K[ and II V Cl without questioning deeply its range 
of validity. However, the smallness of the perturbation is not a sufficient criterion. Indeed, as reviewed in Ref. [23| . 
perturbation theory that relies on linearization has to be questioned, as it yields a solution of the linearized set of 
equations only. Whether it also generates a solution of the full set has to be shown explicitly, and this may be a hard 
task. For instance the full proof of existence of Landau damping [^J] in a Vlasovian frame was a mathematical tour 
de force, the equivalent of a Kolmogorov-Arnold-Moser theorem for continuous systems, and led one of its authors to 
be awarded the 2010 Fields medal. 

As to collisional transport, as usual, it was assumed that some kind of decorrelation justifies the calculation of the 
transport coefficients over the finite times considered in section IVTl However, in reality particle dynamics is chaotic, 
and one is facing the calculation of transport coefficients for a chaotic motion. One should be cautious as, for the 
motion of a charged particle in a spectrum of longitudinal waves, a perturbation calculation yields the quasilinear 
estimate for the diffusion coefficient, while a super-quasilinear regime, a synergetic effect in chaos, is found to exist 
in this chaotic dynamics for intermediate resonance overlap 25]. There, diffusion becomes quasilinear for strong 
resonance overlap, but not because the perturbation calculation becomes valid again ([26J], section 6.8.2 of 12]). 

Furthermore, we used only a very specific part of the fundamental equation ([15^ : the one involving linearization 
and smoothing. It would be interesting to study the effect of the coupling of Fourier components with both co- 
herent and incoherent effects, in particular, to perform the analysis of section |V] by substituting fc^t/(m, cj) with 
— ^ ^^gg exp(— ikm • rjo)Rj{u! — i^mj)- The question arises: is it possible to recover the hole solutions propagating 
near thermal velocity or slower, which are smooth and nonlinear structures satisfying the full nonlinear Vlasov-Poisson 
system (Ref. [l^ and references therein)? 
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